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ABSTRACT 

We review results on rare and forbidden decays of D°, D + , and Df mesons from 
experiments at FNAL. The decay modes studied have two leptons in the final 
state and, if observed, would constitute evidence for flavor-changing neutral- 
current, lepton-flavor-violating, or lepton-number-violating processes. To date, 
no evidence for these decays has been observed and upper limits are obtained 
for their branching fractions. These limits can constrain various extensions to 
the Standard Model. We present new upper limits from FNAL E791 on the 
branching fractions for more than two dozen three- and four-body decay modes. 



1 Introduction 



Searches for rare and forbidden decays of charm are concerned with final 
states containing two charged leptons. Such processes occur via flavor-changing 



neutral-current amplitudes, lepton-flavor-violating amplitudes (leptons belong- 
ing to different families), or lepton-number- violating amplitudes (leptons be- 
longing to one family but having the same sign charge). Diagrams for these am- 
plitudes typically contain new types of particles having high masses; thus, these 
decays probe energy scales which cannot be accessed directly. For example, 
the amplitude for the flavor-changing neutral-current decay D + — ► tt + fi + /j~ 
is expected to be proportional to g 2 /M x x (phase space), where g is a cou- 
pling constant and M x is the mass of some unknown propagator (see Fig. Q). 
The amplitude for the Standard Model decay D + — ► K n + v is proportional to 
g 2 /M$ v x (phase space), and thus r(7r+ fj+ n~) /T(K° \i + v) = g 4 M^/(g A M^) x 
(phase space ratio). If g « g, then 



M x S3 M w 



B(D+ -> K®n+v) 
—, — , : — : r x (phase space ratio) 

B(D+ — > TT+fL+^j 



1/4 



(1) 



Inserting numbers one finds that a branching fraction B(D° — > 7r + ^i + /i ) w 10 

x 



corresponds to a mass M x w 700 GeV/c 2 . 




Figure 1: Nonstandard flavor- changing neutral- current decay (left) and Stan- 
dard Model charged- current decay (right). 

Over the past year, FNAL E791 has published new limits on branching 
fractions for 24 different two- and three-body decay modes @) , and submitted for 
publication new limits for 27 additional three- and four-body decay modes @) . 
In most cases the E791 limits are the most stringent available. Here we review 
these results and also briefly discuss competitive results from other FNAL 
experiments. 



2 The E791 Experiment 



FNAL E79lj] is a hadroproduction experiment studying the weak decays of 
charm mesons and baryons. The charm particles were produced by imping- 
ing a 500 GeV/c ir~ beam on five thin target foils. The most upstream foil 
consisted of platinum; the other foils consisted of carbon (diamond) . All foils 
were separated by about 15 mm such that D mesons decayed predominately 
in the air gaps between foils. The experimental apparatus H) consisted of a 
silicon vertex detector followed by a two-magnet spectrometer, two segmented 
Cerenkov counters for hadron identification, an electromagnetic calorimeter for 
electron identification, and iron shielding followed by scintillator counters for 
muon identification. The downstream silicon vertex detector consisted of 17 
planes of silicon and was used to reconstruct decay vertices downstream of the 
interaction vertex. The spectrometer consisted of 35 planes of drift chambers 
and two proportional wire chambers. The two dipole magnets bent particles in 
the horizontal plane and had p T kicks of +210 GeV/c and +320 GeV/c. The 
Cerenkov counters contained gases with different indices of refraction; together 
they provided ir/K/p discrimination over the momentum range 6-60 GeV/c. 
Data were recorded using a loose transverse energy trigger. After reconstruc- 
tion, events with evidence of well-separated interaction and decay vertices were 
retained for further analysis. The experiment took data from September, 1991 
to January 1992, recording the world's largest sample of charm decays at that 
time. The final number of reconstructed decays is over 200 000. 

3 Event Selection 

E791 has searched for two-, three-, and four-body D° decays such as D° — > 
e + e+ D° — * <j> /i + fi~ , and D° — * K~ir + /i~e + ; and three-body D + and D+ 
decays such as D + — > ir + fi + fi~ and — > K + fi + fi~. Here and throughout 
this paper, charge-conjugate modes are included unless otherwise noted. The 
sensitivity of a search is determined (or normalized) by counting events in 
a topologically-similar hadronic decay channel such as D + — > K~ir + 7r + for 

x The collaboration consists of: CBPF (Brazil), Tel Aviv, CINVESTAV 
(Mexico), Puebla (Mexico), U. C. Santa Cruz, Cincinnati, Fermilab, Illinois 
Institute of Technology, Kansas State, Massachusetts, Mississippi, Princeton, 
South Carolina, Stanford, Tufts, Wisconsin, and Yale. 



Table 1: E791 selection criteria based on tracking and vertexing. 



Selection criteria 


Value 






(z, — Z. 4- (7 2 


>20/12/12 (D + /D + /D°) 


1 (J.UL. LeLIgCL Ui_lglJ 1 ' a^L, 


> 5 


^dec ^last target 


< 16 mm 


Xtrack 


< 5 


Xvertex 


< 6 


D impact parameter (i.p.) 




w/r/t int. vertex 


< 30/40 A*m (L>° 3-body, 4-body/others) 


-p-|- f i.p. w/r/t decay vertex ^ 
lltracks I i.p. w/r/t int. vertex J 


< 0.01 /0.001/0.0005 (2-body/3-body/4-body) 


p T (transverse to D direction) 


< 0.20/0.25/0.30 GeV/c {D+/D+/D ) 


t = m D x (z dec - z int )/p 


< 5/3/(3 or 2.5) ps (D+/D+/D ) 



which the branching fraction is known. For all searches, the event selection 
proceeded via a "blind analysis" technique in order to avoid biasing the choice 
of selection criteria. This technique has three steps: (a) all events having a 
reconstructed mass within a mass window or "box" around m D are removed 
from the sample; (b) the selection criteria are chosen by optimizing the ratio 
S/ y/B, where S is the number of signal events from a Monte Carlo simulation 
that pass all criteria, and B is the number of events from data that are within 
a "background box" which is near -but exclusive of- the signal box; (c) the 
finalized selection criteria are applied to the events within the signal box to 
see if any candidate events remain. The selection criteria resulting from this 
procedure are listed in Tables |l| and |^. The most important criterion is that of 
SDZ, which is defined as the distance between the interaction vertex and the 
decay vertex divided by the error in this quantity. Values used for this criterion 
were 12 for D° and decays and 20 for longer-lived D + decays. 

The mass windows chosen for the various searches depend upon the num- 
ber of particles in the final state and on whether any of them are electrons. 
For final states not containing electrons, the windows extended ±35, ±30, or 
±20 MeV/c 2 around m D for D°, D + , and Df decays, respectively. For final 
states containing one or more electrons, the mass windows were asymmetric, 



Table 2: E791 selection criteria based on particle identification. 



Muon: 



Electron: EMPROB > 90 (calorimeter response consistent with e) 

(1) P^ > 8 GcV/c 

(2) y-paddle within 1 a of track projection is hit 

(3) ^-position as calculated from y-paddle TDC 
time is within 60 cm of track projection 

(4) If y-paddle within 1 a of track projection is not hit, 
then x-paddlc to which track projects is hit 

Pf, , > 0.10 (£>+ -> Knw norm, for £>+ -> ir££ ) 

Coronkov (s) ' 
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Kit-kit norm, for D° — > 7T7t££, 
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■ K££, 

► 07T norm, for Df — » i-Ttt, 
Kiririr norm, for D° — > Kir££ ) 
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extending 40-70 MeV/c 2 farther below m D than above it as the mass distri- 
butions had low-energy tails resulting from final-state bremsstrahlung. For 
search modes containing a p , K , or <fi meson in the final state, it was re- 
quired that |m 7r+7r _ — m | < 150 MeV/c 2 , \m K _ n+ — m K ,\ < 55 MeV/c 2 , or 
\m K+K _ — < 10 MeV/c 2 , respectively. 

After all selection criteria were applied, no significant excess of events 
above the estimated background was observed. The experiment thus sets upper 
limits as follows: 

UL(D X) = f^-) ( £ -f^) x B norm , (2) 

\ iv norm/ V b X J 

where N x is the 90% C.L. upper limit on the mean number of signal events as 
determined from the number of candidate events observed and the estimated 
background; iV norm is the number of events observed (after background sub- 
traction) in a hadronic normalization channel such as D + — ► K~tt + tt + : s x and 
e norm are the overall detection efficiencies for the search channel (D —* X) and 



Table 3: Hadronic decay modes used to normalize the E791 searches for 
D° -► £+£2, D+ s) -► Pt x £ % , D° -► V7+€T, and L>° -► h x h 2 £ x £ 2 decays 
(P = pseudoscalar, V = vector). 



Search mode 


Normalization mode 


AT norm 


D°- 


■+ ^ £ 2 


D°- 




25210±179 


D+ - 


^P£,£ 2 


£>+ - 


-* K~TT + 7T + 


24010 ±166 


Df- 


->Piit 2 




-* (j)TT + 


782 ± 30 


D°- 


♦ pH\£ 2 






2049 ±53 


D°- 


♦ K*°£i£ 2 




-> /V 7T 7T 


5451 ±72 


D°- 








113 ±19 


D°- 


-> irir£ 1 £ 2 




-> 7T + 7T~7r + 7r~ 


2049 ±53 


D°- 


-> if7r^ 2 






11550 ±113 


D°- 


-> KK£ X £ 2 






406 ±41 



normalization channel, respectively; and B nmm is the branching fraction for the 
normalization channel as taken from the Particle Data Book@) . The upper lim- 
its N x are calculated using the method of Feldman and Cousins H) in order 
to account for estimated background. They are subsequently increased via the 
prescription of Cousins and Highland & to account for systematic errors. 

There were eight hadronic decay channels used for normalization. These 
channels are listed in Table |^ along with the number of events obtained for 
each after background subtraction. In general, a dilepton search mode was 
normalized to a Cabibbo-favored hadronic mode having the same number of 
tracks in the final state and, whenever possible, the same daughter particles 
except for the substitution of pions for leptons. For example, the D + — ► K£ x £ 2 
searches were normalized to D + — > K~ir + ir + decays. However, the D° — ► 
p°£f£ 2 searches were normalized to D° — -> ir + tt~ tt + tt~ decays, as there is no 
published branching fraction for D° — > p°w + ir~ . 

4 Background Estimate 

There were two main sources of background in E791: "reflection" background 
arising from fully-reconstructed hadronic D decays in which two of the tracks 
were misidentified as leptons, and "combinatoric" background arising from ac- 



Table 4: Hadronic decay channels contributing reflection background to the 
E791 D° -► l+i 2 , D+ s) -> Pl x l 27 D° , and D° -*• h x h 2 t x t 2 samples 

(P = pseudoscalar, V = vector). 



Cabibbo-favored 


Cabibbo-suppressed 


D° -> K-TT+ 


L>° -> 7T+7T- 


D+ -> K-7T+7T+ 
73+ — > 7T~7T + 7r + 

£>+ -> K-K+TT+ 
A+ — > pK^n + 


Z}+ — > 7r _ 7r + 7r + 
L>+ K-K+TT+ 
D+ K+n+TT- 


D° -> i^-7r+7r-Tr+ 


D° — > 7r + 7r~7r + 7r~ 
£>° -> K+K-TT+n- 



cidental combinations of tracks and vertices. Most of the reflection background 
was eliminated by excluding events with invariant masses (assuming all daugh- 
ters to be 7r or K) near m D ; i.e., it was required that \m(h 1 h 2 h 3 h 4 ) — m D0 \ > 
35 MeV/c 2 , {mfahzha) - m D+ \ > 30 MeV/c 2 , and ^{h^h^} - m fl +| > 
20 MeV/c 2 . These requirements were imposed for all h 1 . . .h 4 final states listed 
in Table ||, with one exception: those final states having the same number of 
kaons as that of the search mode were not excluded in this manner, as the ac- 
ceptance loss for signal events would have been excessive. Instead, background 
from these modes was estimated as follows. 

First, the probability for a pair of pions to be misidcntified as /x/x, fie, 
or ee was estimated from data. This probability was then multiplied by the 
number of events observed in the hadronic channel for which the reflection cut 
could not be applied. A factor was included to account for the fraction of 
these events that would reflect into the signal mass window if two daughter 
pions were misidentified as leptons. For example, the reflection background in 
the D+ — > Tr+fi+fi- sample arising from D+ — » 7r~7r + 7r + decays is calculated 
as P-.-...., x N _ , . x fx J, where f is the fraction of D+ — ► tt~-k+ , k+ 
decays that would reflect into the D+ mass window if two daughter pions were 
misidentified as muons (75% in this example). The factor J c accounts for 
the different ways a hadronic mode can be misidentified as a dilepton mode; 
e.g., there are two ways that D+ — » 7r~7r + 7r + can be misidentified as D+ — > 

7T+ ~[X+ fl~ . 



The misidentification probabilities for three-body decays were obtained 
from the final D + — > K~&ji^ samples, and those for four-body decays from the 
final D° — > K~ir + £ samples, where all candidates observed (after subtract- 
ing combinatoric background estimated from mass sidebands) were assumed to 
originate from D + — > K^tt + tt + and D° — ► K~Tr + ir~Tr + decays, respectively. 
For example, there were 13, 5.2, and 6 events passing all selection criteria for the 
K~p + p + , K~ n + e + , and K~e+e + samples, respectively, and 17 730 events (af- 
ter background subtraction) in the D + — K^tt + t: + sample. Thus, the misiden- 
tification probabilities for three-body decays were P — (7.3 ± 2.0) x 10~ 4 , 
P Me = (2.9 ± 1.3) x 10~ 4 , and P ee = (3.4 ± 1.4) x 10~ 4 . The reflection back- 
grounds for the D + — » K~t[G± and D° — » K~n + €[ £ 2 modes themselves were 
taken to be zero as there was no independent estimate of the misidentifica- 
tion rates. This results in conservative upper limits. Unfortunately, there were 
too many events observed for the D + — > K~ t[£ 2 modes (after accounting for 
combinatoric background - see below) to set meaningful upper limits. 

After the mass reflection cuts, combinatoric background was estimated by 
averaging the number of events in the mass sidebands both above and below 
the D signal mass window and scaling this number by the size of the signal 
mass window relative to that of the mass sidebands. If there were no events 
in the higher mass sideband, it was assumed that there were no combinatoric 
background events in the signal box. This assumption avoids overestimating 
background due to statistical fluctuations. Because it tends to underestimate 
background, it results in a conservative upper limit. 

There were also small backgrounds to the D° — > K~tt + £ + £~ , D° — > 
K+K-£+£-, D° -> K*°£+£-, and D° -> (f>£ + £- samples arising from D° -> 
K-jt+p , D Q -> K+K-p Q , D° -> K*°p°, and D° -> 0p° decays, respectively, 
where p° — > £ + £~. These backgrounds were estimated using the branching 
fractions for D° — > h 1 h 2 p°, D° —> Vp°, and p a — > £ + £~ from the Particle Data 
Book®. 

5 Upper Limits 

The final event samples after all selection criteria were applied and signal boxes 
"opened" are shown in Figs. The number of background events estimated, 
the number of candidate events observed, the overall systematic error, and the 
resultant 90% C.L. upper limits are listed in Tables || and ^. The systematic 



errors arise mainly from four sources: errors resulting from the fits to the nor- 
malization channels, statistical errors on the number of Monte Carlo events 
generated and accepted (used to calculate acceptance), uncertainties in the 
amounts of reflection and combinatoric background, and uncertainties in the 
relative detection efficiencies between the search modes and their normaliza- 
tion channels. The upper limits on the branching fractions are compared to 
(previous) limits from the Particle Data Book E|) in Fig. [|. Of the 51 decay 
modes listed, all but six have upper limits more stringent than previously pub- 
lished results. For 26 of these modes, the E791 limits are the first such limits 
reported. 

6 Other Experiments 

There are five published limits that remain superior to those obtained by E791, 
and one published limit that is equivalent. These were obtained by the CLEO 
experiment 0) (D° -> p°e + e-, p°n ± e T , <j>e + e', ^M ± e T ), by the BEATRICE 
experiment I) (D° -► M + AO> and b Y FNAL E687 W (D+ -► K+e+ e -). This 
last experiment used a photon beam of mean energy ~ 220 GeV to photopro- 
duce charm, and a silicon strip vertex detector (like E791) to reconstruct D 
decay vertices. The experiment ran concurrently with E791 and obtained an 
upper limit for B(D + — > K + e + e~) identical to that from E791. 

Another hadroproduction experiment, FNAL E771 0), also ran concur- 
rently with E791 but used an 800 GeV/c proton beam. This experiment ob- 
tained an upper limit for B(D° — > /i + /i~) remarkably close to that obtained by 
BEATRICE. In fact, the 90% C.L. upper limits from BEATRICE, E771, and 
E791 are within about 25% of each other: B(D° -> fi + ^) < 4.1, 4.2, and 5.2 
xlO -6 , respectively. Assuming these results uncorrelated, we combine them to 
obtain B(D° -> n + fi~) < 1.5 x 10~ 6 at 90% C.L. This limit is still many orders 
of magnitude larger than the Standard Model expectation 0) (dominated by 
long-distance effects) of ~ 10~ 15 . 

For the future, we expect greater sensitivity than that of the above exper- 
iments from FNAL E831 (FOCUS) HI), a photoproduction experiment that is 
an upgraded version of E687. The E831 detector employed more muon count- 
ers than did E791 and also used a finer-grained electromagnetic calorimeter. 
The experiment took data in 1996-97, recording a charm sample approximately 
four times larger than that of E791. The analysis of this data set is underway. 



7 Summary 



In summary, E791 has completed an extensive search for flavor-changing neutral- 
current, lepton-flavor- violating, and lepton- number- violating processes and sees 
no evidence for these decays. The experiment has set upper limits on 51 differ- 
ent decay modes; all limits but six are improvements over previously published 
results. Many of these limits are at the 1CP 5 to 10~ 6 level and can constrain 
various extensions 111) to the Standard Model. We anticipate even more strin- 
gent limits (or possibly signals) from E831, which has a "cleaner" data sample 
due to the photoproduction process, superior muon and electron identification, 
and approximately four times as many reconstructed charm decays. 

References 

1. E. M. Aitala et al, Phys. Rev. Lett. 76, 364 (1996); Phys. Lett. B462, 401 
(1999). 



2. E. M. Aitala et al, FERMILAB-Pub-00/280-E (2000), |hep-ex/0011077 
(submitted to Phys. Rev. Lett.). 

3. E. M. Aitala et al, Eur. Phys. J. direct C4, 1 (1999). 

See also: E. M. Aitala et al, Phys. Lett. B371, 157 (1996); Phys. Rev. 
D57, 13 (1998). 

4. D. E. Groom et al. (Particle Data Group), Eur. Phys. J. C15, 1 (2000); 
C. Caso et al. (Particle Data Group), Eur. Phys. J. C3, 1 (1998). 

5. G. J. Feldman and R. D. Cousins, Phys. Rev. D57, 3873 (1998). 

6. R. D. Cousins and V. L. Highland, Nucl. Instrum. Meth. A320, 331 (1992). 

7. A. Freyberger et al, Phys. Rev. Lett. 76, 3065 (1996); Phys. Rev. Lett. 
77, 2147 (1996). 

8. M. Adamovich et al, Phys. Lett. B408, 469 (1997). 

9. P. L. Frabetti et al, Phys. Lett. B398, 239 (1997). 

10. T. Alexopoulos et al, Phys. Rev. Lett. 77, 2380 (1996). 

11. S. Pakvasa, UH-511-871-97 (1997), |hep-ph/9705397j 



12. P. D. Sheldon, Charm Mixing and Rare Decays, in: Heavy Flavors 8, Proc. 
of the 8th Int. Symp. on Heavy Flavour Physics, Southampton, UK, 25-29 
July 1999 (The Journal of High Energy Physics Conf. Proceedings, eds. 
P. Dauncey and C. Sachrajda, |http : / /jhep . sissa. it ). 

See also: J. M. Link et al, Phys. Lett. B485, 62 (2000); Phys. Lett. B491, 
232 (2000). 

13. See for example: G. Lopez Castro, R. Martinez, and J. H. Muhoz, Phys. 
Rev. D58, 033003 (1998). 



m(7Y + e + e ) 



> 

CD o 



CD 

CL 



D 



D 

CJ 



m(TT 



m(n e + e + ) 



1.8 2 
m(K + e + e~) 



m(TY + /^ + /^ ) 



m(7T + e + e ) 



m(7Y y a + /x + ) 



1 .8 2 
m(7T~e + e + ) 



m(K + /i + /i, ) 



1.8 2 
m(K"/2 + /i + ) 



FOl , Q, 



m(K e + e + ) 




1.8 2 
m(7T _ /x + e + ) 



1 .8 2 

m(K + ,u + e~) 



.in.- -i i -' i-j 



1 .8 2 
m(7T + /i, + e~) 



m(7T a e ) 



1.8 2 
m(K + e + e~) 



1 .8 2 
m(K + /2 + 



1 .8 2 
m(K"/i + e + ) 



1.8 2 1.8 2 

m(,u*,tO m(e*e~) 

Invariant Mass (GeV/c 2 ) 



1.8 2 



Figure 2: final L>+ -> P^ 2 (Vows L>+ -> P^ 2 (Vows ^-7j, and D° -» 

■^i"^2 ( /row ^ ewen ^ samples. The solid curves represent estimated background; 
the dotted curves represent signal shape for an event yield equal to the 90% 
CL. upper limit; the dashed vertical lines denote the signal mass windows. 



> 

2 
m 



























r 




I 


_ 




In 


1 


~n, 




j n 


[1 




1.8 

mD° 




2 




1 .3 

mD° -> 


7T + TT"e + e~ 


2 




1 .8 

mD° 




2 






k\ n ,h n nn 


I 

_ I 

T] ,L JWoj] , 


1 ! fi 
1 n, 1 n wh 


'X™ 




1.8 
mD° 




2 




1 .8 

mD° -> 


K 7T + e + e 


2 




1 .8 
IT D° 




2 


i 
i 

ir . on ! ! u w 


- li ..m < 


I 
I 


_ 


l&A 


n ir 




1.8 
mD° 


-> K~K»~ 


2 




1.8 

mD° -> 


K"K + e + e" 


2 




1 .8 
mD° 




2 






.-- I 


I 


_ 

.j.J ] r 


1 
1 
1 
1 

-.. 1 , 1 


i 

_ i 




1 .8 
raD 


-» pmV 


2 




1.8 

mD° — 


► pe + e 


2 




' .8 
m 




2 


i 

" , i M-. ! w , 


I 
I 

" .! u ! . n 








1.8 

mD° 


— > K'fi + fi 


2 




1.8 

mD° -> 


K"e + e 


2 




1 .8 

mD° — > Kfi*e~ 


2 


1 
1 

" , i ...(•'"'•■._. i n , 


i 
i 

j. j i 


i 
i 
i 

i i i 




1.8 
mD c 


— > <Pf^^ 


2 




1.8 

mD° — 


- (fie + e 


2 




1 .8 

mD° H» <p^e" 


2 


i 
i 

, w i !.--.. n! , n 


i 
i 

_ i 

j.j i--... i 


i 
i 

" n la-, n , 




1.8 
mD° 


TT TT /U + /j- + 


2 




1.8 

mD° -> 


ix 7i e + e + 


2 




1 .8 
mD° 


— > 7T 7T U + e + 


2 


_ 


; \ 

III n,i 


i 
i 

y -,-fl, ! , n 


_ 

k j.!.l..J 






1.8 
mD° 


"> K"7T» + 


2 




1.8 

mD° -> 




2 




' .6 
mD° 




2 


_ 

r , n n n-\ ! n n r 


i 
i 

:...j1... ii ik ! , h.. 


_ 

... .,..1 J , L 




1.8 
mD° 


-> K"K>V 


2 




1.8 

mD° — > 


K"K"e + e + 


2 




1 .8 
mD° 


-> K~K>V 


2 



nvariant Mass (GeV/c 2 ) 

Figure 3: Fmai £>° — ► Vt[l 2 and D° — » h 1 h 2 £ 1 £ 2 event samples for nonreso- 
nant modes (rows 1-3), resonant modes (rows 4~6), and same-signed dilepton 
modes (rows 7-9). The solid curves represent estimated background; the dotted 
curves represent signal shape for an event yield equal to the 90% C.L. upper 
limit; the dashed vertical lines denote the signal mass windows. 



Table 5: The number of background events estimated, the number of candidate 
events observed, the overall systematic error, and the 90% C.L. upper limits on 
N x and the branching fraction for D° — > if £2 an d Df. — > P decays. 
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Table 6: The number of background events estimated, the number of candidate 
events observed, the overall systematic error, and the 90% C.L. upper limits on 
N x and the branching fraction for D° — ► Vlfl^ and D° — > h 1 h 2 £ 1 £ 2 decays. 
The total backround estimate (N tot ) includes small contributions from D° — ► 
K-n+p a , D° -» K+K-p°, D° -» 7?*V, and L>° -» <^>°, w/iere p° 
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